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ABSTRACT 

We present new VLT spectroscopic observations of the most distant quasar known, SDSS J1030+0524 
at z=6.28, which was recently discovered by the Sloan Digital Sky Survey. We confirm the presence of 
a complete Gunn-Peterson trough caused by neutral hydrogen in the intergalactic medium. There is no 
detectable flux over the wavelength range from 8450 to 8710 A. We set a stronger limit on the drop of the 
flux level blueward of the Lyct line: a factor of > 200. Below 8450 A the spectrum shows a rise in flux, 
with a large fraction (> 60 %) of the total emission produced by few narrow features of transmitted flux. 
We discuss the proximity effect around this quasar, with the presence of transmitted flux with many 
absorption features in a region of about 23h~ 1 1 2 3 4 comoving Mpc. If assuming the surrounding medium 
were completely neutral, the size of this region would imply a quasar lifetime of ~ 1.3 xlO 7 years. 

We also present near-IR spectroscopy of both SDSS J1030+0524 and of SDSS J1306+05, the second 
most distant quasar known, at redshift 6.0. We combine measurements of the CIV line and limits on 
the Hell emission from the near-IR spectra with the NV line measurements from the optical spectra to 
derive the metal abundances of these early quasar environments. The results are indistinguishable from 
those of lower redshift quasars and indicate little or no evolution in the metal abundances from z~ 6 to 
z~2. The line ratios suggest supersolar metallicities, implying that the first stars around the quasars 
must have formed at least a few hundreds of Myrs prior to the observation, i.e., at redshifts higher than 
8 . 

Subject headings: cosmology: observations — galaxies: formation — galaxies: quasars: absorption 
lines — galaxies: quasars: individual (SDSS 1030+0524, SDSS 1306+0356) 


1. INTRODUCTION 

Quasars are amongst the most luminous objects in the 
Universe, allowing us to study them and any intervening 
material out to very large distances, corresponding to look- 
back times when the Universe was very young. Hence 
finding and studying quasars at high redshifts is one of 
the best ways to constrain the physical conditions in the 
early Universe. The mere existence of luminous quasars at 
such early times, and the implied presence of black holes 
with M > 10 9 M 0 place stringent limits on the epoch at 
which massive condensed structures formed, thereby con¬ 
straining structure formation models (e.g. Efstathiou and 
Rees 1988). At high redshifts such luminous quasars must 
be associated with very massive halos, hence they are ex¬ 
pected to be found near high peaks (4-5 a or more) of 
the large scale density field (e.g. Efstathiou & Rees 1988, 
Nusser & Silk 1993) that collapsed sufficiently early. Fur¬ 
thermore, several arguments suggest that the formation of 
elliptical galaxies and massive bulges are paralleled by an 
early quasar phase (e.g. Kauffmann & Haehnelt 2000). 


The spectra of high redshift quasars contain important 
information on the enrichment history of the gas in the 
quasar environment, and probe the star formation preced¬ 
ing the epoch at which the quasars are observed, possi¬ 
bly the first stars that formed in massive collapsed struc¬ 
tures (e.g. Hamann & Ferland 1999). Finally, high redshift 
quasars serve as probes of the intergalactic medium via the 
absorption of the Lyman a forest. 

The Sloan Digital Sky Survey (SDSS - York et al. 2000) 
has, amongst its scientific aims, the construction of the 
largest sample of quasars ever, with more than 10 5 ob¬ 
jects spanning a large range of redshift and luminosities. 
The SDSS has already found an unprecedented number of 
new high redshift quasars, including more than 200 new 
quasars at z>4 (e.g. Fan et al. 2000, Zheng et al. 2000, An¬ 
derson et al. 2001, Fan et al. 2001a, Schneider et al. 2001). 
These high redshift quasars have been efficiently selected 
by their distinctive position in color-color diagrams, with 
characteristic colors due to the main feature of the quasar 
spectra, viz., the strong Lya emission line, the Lya forest 
and the Lyman limit. 
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Recently, Fan et al. (2001b) presented the discovery of 
SDSS J1030+0524, found during follow-up spectroscopy 
of i-band drop-out objects, i.e. objects showing very red 
i* — z* color and relatively blue z*-3 colors (following the 
notation of previous SDSS papers, we use the superscript 
* for the photometry, and the letters alone for the filters, 
Stoughton et al. 2002). Moderate signal to noise ratio 
(S/N) optical spectra taken with the ARC 3.5m telescope 
showed very strong (rest-frame equivalent width EW~50 
A) Lya + NV emission at ~8850 A and ~9050 A, respec¬ 
tively. Based on a fit of the CIV line the redshift was esti¬ 
mated estimated to be z~ 6.28 which makes it unambigu¬ 
ously the most distant quasar known. Follow up optical 
spectroscopic observations with Keck (Becker et al. 2001) 
have shown the first clear detection of a complete Gunn- 
Peterson trough (Shklovsky 1964, Scheuer 1965, Gunn & 
Peterson 1965) in the spectrum of this quasar. The flux 
level drops by a factor of at least 140 relative to an es¬ 
timate of the unabsorbed continuum level, and the spec¬ 
trum is consistent with zero flux in the Lya forest region 
immediately blueward of the Lya emission line. Even if 
the existence of the trough itself does not imply that the 
quasar is observed prior to the re-ionization epoch (e.g. 
Barkana 2001), the fast evolution of the mean absorption 
in the spectra of quasars at z> 5 suggests that the universe 
is approaching the reionization epoch at z~ 6 (Becker et 
al. 2001, Fan et al. 2001b, Djorgovski et al. 2001). In a 
companion paper (Fan et al. 2001c), we use cosmologi¬ 
cal simulations to estimate the evolution of the ionizing 
background and to constrain the redshift of reionization 
(see also Gnedin 2000, Cen & McDonald 2001, Lidz et al. 
2001 ). 

We have obtained an optical VLT spectrum of the 
quasar to independently measure the Gunn-Peterson 
trough and to set tighter limits on the transmitted flux. 
We have also obtained a near-IR spectrum, to set con¬ 
straints on the metallicity of the quasar environment and 
to better estimate the redshift from the CIV emission line, 
which unlike Lya is not affected by intervening absorption. 
Finally we obtained a near-IR VLT spectrum of the second 
most distant quasar SDSS .11306+0356 at z=5.99 (Fan et 
al. 2001b). We will use the results from this spectrum in 
the discussion of the metallicity in high redshift quasars. 

The paper is organized as follows: in Section 2 we de¬ 
scribe the observations and present the optical and near- 
IR spectra of the two quasars. In Section 3 we discuss 
the main results obtained from the optical spectrum: the 
Gunn-Peterson trough and the limits we can set from the 
new data, as well as the region in which flux is present, and 
the proximity effect. In Section 4 we derive the metallicity 
in the two highest redshift quasars and discuss its impli¬ 
cation for the star formation timescales. 

Throughout the paper we assume a flat, A-dominated uni¬ 
verse with Hq = 65 km s -1 Mpc -1 , A = 0.65 and 0 = 0.35 
(Ostriker & Steinhardt 1995; Krauss & Turner 1995). 

2. OBSERVATIONS AND RESULTS 
2.1. Optical spectroscopy 

Optical long slit spectroscopic observations were carried 
out in service mode on 28 and 30th of June 2001, with 
the FORS2 spectrograph mounted on VLT/YEPUN on 
Paranal. On each night, the quasar SDSS 1030+0524 was 


observed for a total of 3600 sec divided into two 30 minute 
exposures. The target was placed at different locations 
along the slit for each individual frame, to minimize prob¬ 
lems related to flat fielding and other possible systematic 
errors. Weather conditions were clear but not photomet¬ 
ric. The seeing was of order 0.8" on both nights. We 
used a 1" slit together with the new holographic grism 
GRISM_1028z+29 and the order separation filter OG590, 
which has greatest sensitivity at A « 8800A, and covers 
the wavelength range we are interested in (7900 A to 9300 
A) This grism also offers the highest resolving power, R~ 
2700, corresponding to a resolution of ~ 3.1 A, which will 
allow us to cleanly resolve out the night sky, thereby sig¬ 
nificantly increasing our signal-to-noise ratio and our abil¬ 
ity to subtract the sky. The dispersion of the spectrum 
is 0.68 A per pixel. Given the rather high airmass of 
1.4 < sec(z) < 1.9 for the observations, the slit was kept 
at the parallactic angle. 

Image reduction (including bias subtraction, flat fielding 
and correction for slit distortion for each individual frame) 
was performed with standard IRAF procedures. The most 
critical part of the reduction was the sky-background sub¬ 
traction, especially due to the presence of strong sky emis¬ 
sion lines in the wavelength region where the flux is very 
close to zero. In particular the absolute level of flux de¬ 
tection changes slightly (by ~ la) depending on the pa¬ 
rameters used for sky subtraction. We therefore assumed 
that outside the trace of the spectrum of the quasar, i.e. 
at other places along the slit, there was no flux. Therefore 
after the normal sky subtraction, we determined the av¬ 
erage of the residual sky signal in a large region on either 
side of the spectrum and subtracted this single value from 
the entire data frame. 

Wavelength calibration was done using observations of 
a Xe-Ar lamp obtained during the morning calibration: 
the intrinsic accuracy of the wavelength fit is better than 
0.2 A. To this error we have to add the uncertainty due 
to the large range of airmasses of the observations, which 
we estimate from the spectra as less than half a pixel. 
The total uncertainty is then 0.4 A. Relative flux calibra¬ 
tion was performed using a long slit spectrum taken with 
the same setup of the spectrophotometric standard EG274 
(Hamuy et al. 1992). The spectra from each night were 
calibrated separately and then coadded using a cosmic-ray 
rejection algorithm. Finally the spectrum was extracted: 
since we could not trace it from the image itself, due to 
the absence of intrinsic quasar flux throughout the Gunn- 
Peterson trough (see below), we used the trace derived 
from the standard star. The spectrum was not corrected 
for telluric absorption. The two-dimensional spectrum and 
the extracted one-dimensional spectrum (smoothed to 3 
Aper pixel) are presented in Figure 1. The main charac¬ 
teristics of the Lya and NV emission lines are given in 
Table 1. 

2.2. Near-IR Spectroscopy 

Near-IR observations of SDSS 1030+0524 were carried 
out in service mode with ISAAC on VLT-ANTU on 18 
May 2001. Observations were made with the SW Rock¬ 
well 1024x1024 pixel Hawaii array which has a scale of 
0.148 arcsec/pixel and the lower resolution grating, whose 
resolving power with a 0.6" slit is 860. The seeing dur- 
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ing the night varied between 0.4 " and 0.7". The total 
integration time was 5760 sec, split into 2 sets of 16 im¬ 
ages. We used a detector integration time DIT=180 sec 
over which time a single integration was done (NDIT=1). 
The object was moved between exposures by 20" along the 
slit in a ABBA sequence, with an extra jitter of amplitude 
5-10" around each A and B position. The airmass varied 
between 1.15 and 1.27. 

The ISAAC data were reduced in the standard way with 
IRAF, subtracting each image from its adjacent one. Af¬ 
ter flat fielding, we corrected each frame for geometric dis¬ 
tortions using the sky lines as well as a spectrum of an 
argon lamp, which also led to wavelength calibration. Fi¬ 
nally the frames were aligned using the shifts from the 
image headers and combined with the IRAF/avsigclip re¬ 
jection algorithm. The spectrum was then extracted using 
a 3 pixel aperture, with a trace derived from the spec¬ 
trum itself. Unfortunately the telluric star from the ESO 
standard calibration plan was observed with a different 
instrument configuration. Since no telluric spectra were 
present in the ESO database, the observations were re¬ 
peated on the 15th of June with a seeing of around 0.5" 
and at similar airmass as the science observations. The 
standard telluric star (H 78530) was reduced in a simi¬ 
lar way as the spectrum, and in particular the wavelength 
calibration was done with a lamp exposure taken on the 
same night. However, it is apparent that the removal of 
telluric features from our data is not perfect and the re¬ 
sulting spectrum is noisier than we had hoped. Finally the 
spectrum was flux calibrated using the total J band mag¬ 
nitude of the object (J=18.87± 0.10, Fan et al. 2001b) and 
is presented in Figure 2. The measured characteristics of 
the CIV line and the 3er limits on the undetected Hell 
line are reported in Table 1. Contrary to what we hoped, 
the redshift we derive from the CIV line is actually less 
accurate than that derived from the lines in the optical 
spectrum: the reason is the uncertainty in the removal of 
some deep telluric features at the same wavelength as the 
CIV emission line. 

An infrared spectrum of the quasar SDSS J1306+0356 
at z=5.99 was obtained with ISAAC with the aim of de¬ 
tecting the CIV line. The observations were carried out in 
service mode on 1 July 2001; the total integration time was 
5760 sec, which was divided into 2 sets of 16 images, each 
with 180 sec integration time. We used the SZ filter, which 
gives a total wavelength coverage from 0.98 fim to 1.14/zm. 
The airmass during the observations varied from 1.19 to 
1.39 and the typical seeing was 0.8". The observational 
strategy was the same as for SDSS 1030+0524. Reduction 
was carried out as described above. The telluric features 
were eliminated with observation of the telluric standard 
H74113 , observed with the same set-up of the target dur¬ 
ing the same night. The resulting spectrum is shown in 
Figure 3. The two absorption features at A ~9900 A that 
are observed in the spectrum are due to Mgll absorption 
at redshifts 2.20 and 2.53, respectively. No obvious inter- 
galactic CIV absorption is seen. 

From a fit to the peak of the CIV line we derive a 
redshift of 6.00+0.01 for this quasar, consistent with the 
Zem=5.99±0.02 derived by Becker et al. (2001) from the 
optical emission lines (NV). Note that the CIV line is of¬ 
ten offset relative to the systemic quasar redshift, e.g. for 


z~ 2 SDSS quasars Richards et al. (2002) find an average 
shift between the CIV line and the systemic redshift on 
the order of 800 krn/s, with some shifts as large as 3000 
krn/s (see also Vanden Berk et al. 2001, Laor et al. 1995). 
However in our case the CIV redshift agrees within the 
errors with the redshift derived by NV which is normally 
very close to the systemic redshift, so we believe there are 
no such extreme offsets; in the error quoted above we have 
included a possible uncertainty of 500 km/s. The pro¬ 
file of the CIV line actually suggests the presence of both 
broad and narrow components; fitting the overall shape 
with two gaussians gives FWHM for the two components 
of 1600+200 km/s and 4400+400 km/s respectively; the 
broad component is blueshifted relative to the narrow one. 
No limits can be derived for the undetected Hell line since 
it falls at the edge of the spectrum where the S/N is very 
low. 

3. RESULTS 

3.1. The Gunn-Peterson trough: the limits from the VLT 

spectrum 

In the bottom panel of Figure 1, we show the two- 
dimensional spectrum of SDSS 1030+0524 in the wave¬ 
length region between 7990 A and 9250 A. Clearly there 
is a sharp drop in flux, blueward of the Lya emission, and 
the flux is consistent with zero in the region between 8450 
A and 8710 A. 

Following previous papers (e.g. Fan et al. 2001b) we 
have computed the transmitted flux ratio over 3 small red¬ 
shift windows in z abs as T{z abs ) = (f? bs /ff on ), where f° bs 
is the observed flux and ff ont is the intrinsic quasar con¬ 
tinuum flux. For the latter quantity we have assumed a 
power law, with ff on cx f~ 0 - 5 (or f{ on oc fff 1 ' 5 ). We 
normalized this to the continuum flux measured at an ob¬ 
served wavelength of 9200 A (with the normalization value 
f\ = lx 10 -17 erg s _1 cm -2 A -1 ), since this is the only 
range of our optical spectrum free from emission and ab¬ 
sorption lines. 

The average flux ratio in the three wavelength bins is 
0.046+ 0.005, 0.029+ 0.004 and 0.001 ± 0.004 respectively 
in the intervals 8000-8200 A, 8200-8450 A and 8450-8710 
A, approximately corresponding to z abs of 5.65, 5.85 and 
6.05. The error bars include only the photon noise, and do 
not reflect the systematic error due to the unknown contin¬ 
uum shape (which however should be smaller, as argued by 
Fan et al. 2001b). The results are consistent with those re¬ 
ported by Becker et al. (2001) from the Keck spectroscopic 
observations. Therefore, we confirm the non-detection of 
flux found in Keck spectrum in the same wavelength range. 
We emphasize that these are completely independent mea¬ 
surements, using a different telescope at a different site 
with a different instrument, and both spectra show the 
same non-detection of flux. 

The implied drop of flux for the wavelength region 8450- 
8710 A is a factor of 200 at a lcr level, a stronger limit than 
previously reported by Becker et al. (2001), who derived 
a factor of 140, also at lcr. The corresponding effective 
Gunn-Peterson optical depth, r e // = — ln(T), is > 5.3. 
Note that we cannot compute the D^ parameter defined 
in Oke & Korycansky (1982) as D A = (1 - f° bs /ff on ) 
measured in the region between restframe Lya and Ly (3 
(A=1050 and 1170 A) and is normally used to describe the 
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Fig. 1.— Upper panel: the VLT/FORS2 one-dimensional optical spectrum of the quasar SDSS 1030+0524 shown in the observed frame and 
smoothed to a dispersion of 3 A pixel -1 . The bottom line shows the error array derived from the sky spectrum and offset by —0.2 x 10 -17 cgs 
for clarity. Lower panel: a grey-scale representation of the sky-subtracted two-dimensional spectrum plotted on the same wavelength scale. 
The spatial extent of the plot is 40 pixels, corresponding to 20 arcseconds. 
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Fig. 2.— A composite of the optical and near IR (J band) spectra of SDSS 1030+0524. The dashed line is the continuum normalized at 
12400 A. 



Fig. 3.— The z-band spectrum of the quasar SDSS 1306+0356, showing the CIV emission line. 
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decrement of flux below the Lya line, because the wave¬ 
length range covered by our observations is not quite large 
enough. 

The implication for the ionization state of the IGM are 
discussed in Fan et al. (2001c) where it is shown that the 
observed properties of the IGM at z~6 are typical of those 
in the era at the end of the overlap stage of reionization 
when the individual HII regions merge. Thus this obser¬ 
vation suggests that z~6 marks the end of the reionization 
epoch. 

3.2. The spectrum below 8450 A: IGM absolution 
between z=5.58 and 5.95 

The extracted spectrum shows some emerging flux be¬ 
low 8450 A. However, from an inspection of the two- 
dimensional spectrum it appears that a large fraction of 
the flux comes from a few isolated transmission features, 
while the rest of the spectrum remains consistent with zero 
flux. In other words, a simple description of the transmit¬ 
ted flux in this region as arising from a uniform IGM is not 
sufficient. The presence of these gaps of transmitted flux 
can be explained, in the context of reionization, as due to 
cosmological HII regions around the sources that reionize 
the IGM before the epoch of overlap. The characteristics 
of these gaps, such as their sizes and probability distribu¬ 
tion are important since they can put constraints e.g. on 
the nature of the ionizing sources (quasars and/or galaxies) 
on their lifetimes etc. In particular the reader is referred 
to Miralda-Escude et al. (2000) for a detailed discussion of 
gaps in the Gunn-Peterson trough due to individual HII 
regions (see also Loeb & Barkana 2001, Djorgovski et al. 
2001, Haiman & Loeb 1999). 

In Figure 4 (left panel) we show the significance of any 
possible flux detection in the region between 8000 A and 
8750 A: for each resolution element (the spectrum was 
binned to a resolution of 3 A per pixel) the observed flux 
has been divided by its corresponding error array. There 
are no significant features (> 3<r) in the range between 
8400-8750 A, while there are several such features present 
in the shorter wavelength part of the spectrum. All the 
features above 3cr are also clearly recognized in the two- 
dimensional spectrum, and are visible on both spectra 
from the two different nights. We have tried to quantify 
what fraction of the total flux below 8400 A is accounted 
for by these isolated features. To detect features in a sys¬ 
tematic way, we used SEXtractor (Bertin & Arnouts 1996) 
on the two-dimensional sky-subtracted spectrum. We con¬ 
sidered as an emission feature anything detected in more 
than 3 adjacent pixels with a flux larger than ler plus the 
local background as determined in an annulus around the 
feature. We considered only those features that have a 
peak of emission along the trace of the spectrum (± 1 
pixel); furthermore we eliminated all those features that 
were coincident with sky lines (by using the error array of 
the extracted spectrum), since these could be due to resid¬ 
uals of sky subtraction. This is a conservative approach 
since we might eliminate true features. Eight such fea¬ 
tures are detected from SEXtractor. We then measured 
the flux inside the features, assuming the width as given 
by SEXtractor in the dispersion direction, and 3 pixels in 
the spatial direction. The values measured are reported 
in Table 2, together with the central wavelength and the 


width. The flux from these features accounts for 60% of 
the total flux in the 8000-8400 A range. 

To allow comparison with simulations, we also present 
the cumulative distribution of the transmitted flux ratio 
of the region between 8000 and 8450 A (corresponding to 
a mean redshift of absorption z a f, s =5.74) in Figure 4. For 
this plot the spectrum was rebinned to the resolution of 
3 A/pixel, and then divided by the assumed unabsorbed 
continuum flux estimated in the previous section. Ninety- 
nine per cent of the spectrum has a transmitted flux ratio 
less than 20% of the continuum, and the distribution goes 
to 100% very rapidly (note that all results depend on the 
resolution of our spectrum). For comparison, we also draw 
the same distribution for a random realization of a spec¬ 
trum having the same resolution and the same average 
noise as the real data, but with zero average flux. One 
should actually do this pixel-to pixel using the individ¬ 
ual pixel errors (see Fan et al. 2001c) but such a detailed 
analysis is beyond the scope of this paper. In any case the 
“random” line is clearly above the data, implying that the 
observed spectrum has finite non-zero average flux in this 
wavelength region and cannot be solely due to the noise 
alone. McDonald & Miralda-Escude (2001) present the de¬ 
tailed flux distribution from simulated spectra at slightly 
lower absorption redshift. Their prediction could be eas¬ 
ily compared with these data. For example at z a bs =5.2 
they predict that only 3% of the spectrum should have 
a flux higher than 50% of the continuum, and less than 
0.5% should have a flux greater than 75% of the contin¬ 
uum. Simulations show that because the IGM is highly 
clumpy, most of the flux that we observe is transmitted 
from a few of the most underdense voids, consistent with 
what is seen here (see the detailed discussion in Fan et al. 
2001c). 

3.3. The proximity effect 

The spectrum shows excess transmitted flux immedi¬ 
ately blueward of the Lya line, as well as several absorp¬ 
tion lines. This is due to the fact that luminous quasars 
such as SDSS 1030+0524 ionize the surrounding regions 
and create HII regions of radius several Mpc, the so-called 
“proximity” effect (e.g. Carswell et al. 1982). The trans¬ 
mitted flux extends over about 100 A blueward (i.e., to 
8750 A) of the Lya emission peak at 8850 A. This corre¬ 
sponds to a distance of about 23h _1 Mpc comoving, or a 
proper size of D 0 f, s =3.1h _1 Mpc at the epoch z=6.28. 

The size of the HII region associated with the quasar 
depends on the luminosity and lifetime of the quasar and 
the dumpiness of the gas near the quasar. As pointed out 
by Madau & Rees (2000), if the lifetime of the quasar is 
shorter than the recombination time recombination can be 
neglected, and the evolution of the HII region can be de¬ 
coupled from the expansion of the universe. This holds at 
z=6.28 on Mpc scales if the dumpiness factor C is less than 
~10. The photons will travel freely in the highly ionized 
bubble and will be absorbed in a transition layer (I-front). 
This I-front will initially expand at a velocity very close 
to the speed of light, but the Lya and continuum photons 
from the quasar cannot “catch up” with the boundary of 
the HII region until its expansion speed slows down (Cen 
& Haiman 2000), when it first begins to be limited by the 
number of ionizing photons. 
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Fig. 4.— Left: an enlargement of the region between 8000 A and 8750 A, plotted in terms of significance, i.e. at each pixel the observed 
spectrum has been divided by the error array. The dashed line represent the 3 a level: clearly there are no peaks beyond 8400 A. The 
spectrum was binned to the resolution (~ 3 A). Right: the cumulative distribution of the transmitted flux ratio for the spectrum in the 
8000-8400 A region. The spectrum has been rebinned to the resolution ~ 3A/pixel. The dashed line indicates the same distribution for a 
random realization of a spectrum with zero average flux and the same average noise of the real data. 


In practise, the volume of ionized region that we observe 
will be proportional to the number of Lyman-continuum 
photons that are emitted over the source lifetime. In a 
homogeneous medium we can write: 



tQN ph 

n H f 


where N p h is the photo-ionization rate, and nn is the 
mean hydrogen density within Dj. Following the papers 
cited, and in particular Cen & Haiman 2000, we use the 
mean hydrogen density given by tih = 1-6 x 10~ 7 (1 + 
z) 3 {flbh 2 /0.02) cm -3 . However as we remarked earlier, lu¬ 
minous quasars probably reside in significantly overdense 
regions: the unknown bias factor will contribute to the 
uncertainty of the proximity effect region. In the equa¬ 
tion above, f is the fraction of neutral hydrogen present 
when the quasar switches on. From the observations the 
only constraint we have is that / > 0.01 at 2 ~ 6 (Fan 
et al. 2001c). Simulations (e.g. Gnedin 2000) indicate 
that the ionizing fraction could be / = 0.1 or more at 
2 = 6.28. Taking the quasar luminosity as oc i' -a , 
we derive a total photo-ionization rate for the quasar of 
Nph = 2/a x 10 57 s _1 . Therefore the quasar must have 
been on for tg ~ 1.3/a x 10 7 yrs. This is the total time 
during which the quasar has been emitting ionizing pho¬ 
tons, and not necessarily the time elapsed since it first 
turned on. The derived time is of the same order as the 
e-folding timescale for a black hole accreting with 0.1 effi¬ 
ciency of the Eddington luminosity of 4x 10' yr (Salpeter 
1964). 

Note that tg is shorter than the light-crossing time of 
the HII region (~1.5 xlO' yrs). The reason is that photons 
along the line of sight have been emitted at different red- 
shift and therefore at different lookback times. Therefore 
it is not correct to set a minimum lifetime of the quasar 
by saying that those photons that are ionizing the edges of 
the HII region needed a certain time to propagate there, 
travelling at the speed c, since they could have been emit¬ 
ted at a later time. The light-travel argument would only 
be right if we could observe the HII region seen around the 
quasar in the plane of the sky. 

The above calculation does not take into account other 
effects such as the presence of Lyman limit systems near 
the quasar, which would completely absorb the flux and 
produce a cutoff in the spectrum. In fact such a single 
overdense (and therefore neutral) cloud along the line of 


sight between the quasar and us is quite likely to exist close 
enough to the quasar to give a biased answer for the size 
of the ionized region (see Fan et al. 2001c). The presence 
of a Lyman limit system close to the quasar would lead to 
an underestimate of the quasar lifetime. 

The transmitted flux profile blueward of the Lya line 
also shows numerous absorption features which are pro¬ 
duced by the density fluctuations of the medium around 
the quasar but are not resolved in the present spectrum. 
Cen & Haiman (2000) show that much higher spectral res¬ 
olution measurements (R~10 4 ) of the blue side of Lya in 
bright quasars would provide direct estimates of the den¬ 
sity fluctuations of gas on small scales. 

The proximity effect is one of the primary techniques 
used to estimate the photo-ionization rate of the IGM 
(e.g. Bajtlik, Duncan, & Ostriker 1988, Bechtold 1994), 
by determining the distance from the quasar at which the 
optical depth of the Lya forest is half that in the Gunn- 
Peterson trough. However, in the case of our quasar, the 
ionizing background is dominated by the flux of the quasar 
itself, so the profile of the Lya line is insensitive to the ion¬ 
ization state of the IGM. 

4. METALLICITY IN THE TWO HIGHEST REDSHIFT 
QUASARS 

The metallicity ( Z ) of the gas around high-redshift 
quasars contains important information on the enrichment 
history of the gas in the quasar environment, and probes 
the star formation that came before the epoch of the ob¬ 
servations, possibly the first stars that formed in massive 
collapsed structures (e.g. Harnann & Ferland 1999). In 
particular, line ratios involving nitrogen are valuable trac¬ 
ers of chemical enrichment, since this element is selectively 
enhanced by secondary processing in stellar populations 
and its abundance increases roughly as Z 2 . 

There is increasing evidence from measurements of vari¬ 
ous line ratios that the metallicity of the broad line regions 
around quasars have solar or supersolar values even at red- 
shifts higher than 4, implying that there are stars at that 
redshift with ages exceeding 1 Gyr (e.g. Elston et al. 1994, 
Thompson et al. 1999). Clearly, dating the first star for¬ 
mation becomes increasingly interesting at redshifts of 5 
or beyond, as the age of the universe approaches the en¬ 
richment timescale itself. At redshifts 5.99 and 6.28 the 
universe is only 910 Myr and 860 Myr old respectively, in 
our adopted cosmology. 
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Fig. 5.— Measured NV/CIV and NV/Hell line ratios versus redshift for quasars up to z~ 4 reproduced from Hamann and Ferland (1999) 
and with the addition of the new values measured for the 2 highest redshift quasars (represented as triangles). 


From our spectrum of the z=6.28 quasar, we are able 
to measure the characteristics of the Ly a, NV and CIV 
emission lines, which are reported in Table 1. We do not 
detect Hell, but we can set an upper limit based on the 
S/N of the continuum. For the redshift 5.99 quasar, we 
combine our measurement of the CIV line with the Lyct 
and NV values derived by Becker et al. (2001). The derived 
line ratios for the redshift 6.28 quasar are NV/CIV= 0.35, 
slightly lower but consistent with that previously reported 
by Becker et al. (2001), and NV/Hell > 4.3. For the red¬ 
shift 5.99 quasar, the restframe EW of the CIV line is 37± 
6 A; using the reported restframe EW of NV from Becker 
et al. (EW=17.9 A), we derive for this object a NV/CIV 
ratio of 0.67. Note that in all cases we have derived the 
line ratios using the EW of each line and assuming that 
the continuum has the form of a power law with f v oc v~ a 
and a = 0.5 as previously used. 

In Figure 5 we show the measured values for NV/Hell 
and NV/CIV, in a plot of line ratios versus redshift re¬ 
produced from Hamann & Ferland (1999), which includes 
measured values for quasars at different redshifts. The line 
ratios of z~6 quasars appear fully consistent with those 
at redshift between 2 and 4, indicating little evolution in 
the abundances, and hence metallicities, of quasar environ¬ 
ments. However metallicity might be actually more a func¬ 
tion of local density than of redshift (e.g. Cen & Ostriker 
1999); we also know that metallicity scales with quasar 
luminosity, and since the two quasars are the brightest 
objects at that redshift, we might argue that the most lu¬ 
minous objects at 6 have lower metallicity than their 
counterparts at z~ 4. Clearly more statistics are needed 
to determine whether we are starting to see a decline in 
metallicity as compared to z~4. 

Figure 5 also shows the predictions based on the 
Hamann & Ferland (1993,1999) chemical evolution mod¬ 
els. The models that best reproduce the observed ratios for 
the two z> 6 quasars are those of “Giant Elliptical” nitro¬ 
gen enrichment; these models involve short star formation 
timescales (0.5 Gyr) and IMF favoring high mass stars. 
The solar neighborhood models severely under-predict the 
two line ratios (NV/Hell and NV/CIV). If we plot the 


observed line ratios into Figure 6 of Hamann & Ferland 
(1999), we derive a metallicity of roughly 7-8 Z Q for our 
two quasars. We emphasize that the actual value of the 
metallicity is somewhat uncertain, due to the uncertainty 
in deriving reliable abundances from the emission line ra¬ 
tios: e.g. varying the shape of the EUV continuum, vary¬ 
ing the mix of ionization parameter and column density 
in the radiating regions, and allowing for scattering con¬ 
tributions to NV 1240 (e.g. Krolik & Voit, 1998) can all 
perturb the relative line strengths without indicating any¬ 
thing about the abundances. However the main result is 
that the observed ratios definitely require Z > Zq. Note 
that metallicity can be derived self-consistently from the 
two separate line ratios only if Noc Z 2 . 

From the models of Hamann & Ferland (1993) we see 
that the minimum time needed for a stellar population to 
achieve the observed metallicity in any of their models is 
about 300 Myr. In the assumed cosmology the implied 
redshift of star formation is z~8.7, an epoch when the 
universe was about ~560 Myr old. Even if the quasar 
environments had only solar metallicity, the time required 
for any stellar population to reach such a metallicity would 
be on the order of 100 Myrs or more, implying that the 
star formation must have begun well before the observed 
quasar activity. Again these results are effected by the un¬ 
certainties of the models outlined above: other elements 
and in particular Fe, might provide better “clocks” for 
constraining the age of quasars and the epoch of first star 
formation. Fe has a large delayed contribution from type 
la supernovae, with a predicted time delay of about 1 Gyr 
after the initial starburst (Greggio & Renzini 1983, Mat- 
teucci & Greggio 1986) which does not depend on other 
parameters such as star formation rate etc. Therefore the 
ratio of Fe to a elements such as Mg is an absolute clock 
for constraining the ages of star forming regions. This will 
be pursued with follow-up observations in K-band where 
there are interesting features such as the broad FeIIA2960 
and the MgIIA2798 lines. 

5. CONCLUSION 

We have presented new optical and near-IR spectro¬ 
scopic observations of the two most distant quasars known, 
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at redshifts 6.28 and 6.0. We have confirmed the presence 
of a complete GP trough in the highest redshift object, 
derived new limits on the transmitted flux and optical 
depth in this region, and analyzed the characteristics of 
the flux emission in the region. In Fan et al. (2001c) we 
use semi-analytic models to study the evolution of the ion¬ 
izing background and the epoch of reionization from the 
z~ 6 quasars. 

From the optical and near-IR spectra, we have estimated 
the metallicities of these early quasar environments, find¬ 
ing supersolar values for both objects. These high metal¬ 
licities imply that the first stars around the quasars must 
have formed at least a few hundreds of Mpc prior to the 
observation, i.e. at redshifts higher than 8. Follow-up 
observations of other metal emission lines will help us to 
better constrain the ages of the star forming regions. 
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Table 1 

Emission line properties. 


Line 

redshift 

flux 

10~ 15 erg s -1 cm -2 

FWHM 
km s -1 

EW a 

A 

Lya 

- 

2.70 

4000±400 

41.9±4 

NV 

6.29 ±0.01 

0.86 

2400±260 

14.3±4 

CIV 

6.27 ±0.02 

3.2 

5600±400 

42.5±5 

Hell 

- 

<0.20 

— 

< 3.3 


“Restframe EW 


Table 2 

Characteristics of emission features in the 8000-8400 A range of the spectrum. 


mnn 

Ao 

A 

Flux 

10 _17 erg s -1 cnr -2 

Width(lcr) 

A 

i 

8000.4 

1.59 

8.2 

2 

8006.3 

0.43 

2. 

3 

8012.7 

0.92 

4.8 

4 

8047.1 

1.02 

7.5 

5 

8091.4 

1.12 

5.4 

6 

8150.6 

0.49 

4.1 

7 

8253.1 

0.67 

4.8 

8 

8332.9 

1.02 

4.8 



